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The isospin dependence of the spin-orbit potential is inves¬ 
tigated for an effective Skyrme-like energy functional suitable 
for density dependent Hartree-Fock calculations. The mag¬ 
nitude of the isospin dependence is obtained from a fit to 
experimental data on finite spherical nuclei. It is found to 
be close to that of relativistic Hartree models. Consequently, 
the anomalous kink in the isotope shifts of Pb nuclei is well 
reproduced. 

PACS numbers : 21.60Cs, 21.10.Dr, 21.10.Gv 


The Hartree-Fock approach based upon phenomeno¬ 
logical density dependent forces |lj,|],D has proved to be 
very successful in the microscopic description of ground 
state properties of nuclear matter and of finite nuclei 
over the entire periodic table. In all these calculations 
the spin-orbit potential has been assumed to be isospin 
independent. Its only parameter, the strength, is usu¬ 
ally adjusted to the experimental spin-orbit splitting in 
spherical nuclei like 16 0 or 208 Pb. The exchange term, 
however, causes for nuclei with neutron excess a strong 
isospin dependence of the corresponding single-particle 
spin-orbit held. 

In recent years Relativistic Mean Field (RMF) theory 
J|] with nonlinear self-interactions between the mesons 
has gained considerable interest for the investigations of 
low-energy phenomena in nuclear structure. With only a 
few phenomenological parameters such theories are able 
to give a quantitative description of ground state prop¬ 
erties of spherical and deformed nuclei HU at and away 
from the stability line. In addition, excellent agreement 
with experimental data has been found recently also for 
collective excitations such as giant resonances Q| and for 
twin bands in rotating superdeformed nuclei Q . In many 
respects the relativistic mean-held theory is regarded as 
similar to the density-dependent Hartree-Fock theory of 
the Skyrme type PJlc|| 

Recently, however, detailed investigations of high pre¬ 
cision data on nuclear charge radii in Pb isotopes [11 
and of shell effects at the neutron drip line [13 
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have shown considerable differences between the Skyrme 
approach and the relativistic mean held theory. In 
fact, density dependent Hartree-Fock calculations with 
Skyrme |nj or Gogny forces ju|, which used so far 
antisymmetrized, isospin independent spin-orbit interac¬ 
tions, were not able to reproduce the kink in the isotope 
shifts of Pb nuclei (see Fig. [I|). On the other hand, this 
kink is obtained in the RMF theory without any new ad¬ 


justment of parameters | 1 J]. Another considerable differ¬ 
ence has been found in theoretical investigation of shell 
effects in very exotic Zr-isotopes near the neutron drip 
line: in conventional non-relativistic Skyrme calculations 
the shell gap at isotope 122 Zr with the magic neutron con¬ 
figuration N = 82 is totally smeared out 1131, whereas rel¬ 
ativistic calculations using various parameter sets show 
at N=82 a clear kink in the binding energy as a function 
of the neutron number [Q. This difference is caused by 
the different spin-orbit splitting of the single particle lev¬ 
els in these nuclei. Mass calculations within the Finite 
Range Droplet Model (FRDM) , which are based on 
an isospin dependent spin-orbit term carefully adjusted 
to experimental data, are in excellent agreement with the 
relativistic predictions. 

This gives us a hint, that one does not really 
need full relativistic calculations in order to understand 
these differences between conventional density dependent 
Hartree-Fock calculations and RMF theory. A spin-orbit 
term with a properly chosen isospin dependence might 
represent the essential part of a relativistic calculation. 
In this letter, we therefore explore the isospin dependence 
of the spin-orbit term in non-relativistic Skyrme calcula¬ 
tions and analyze the consequences of this on the nuclear 
properties. We start from the Skyrme type force 


F(l, 2) = t 0 (l + x 0 P a )S(r 1 - r 2 ) 

+ ti(l + XiP a )(5(r 1 - r 2 )k 2 + h.c.) 

+ t 2 (l + X 2 P <T )k 6 (ri - r 2 )k ( 1 ) 

+ ^h{l + x 3 P a )5{v 1 -v 2 )p a 

+ W 0 (l + x w P t ){<j^ + ff ( 2 ))k x <5(n - r 2 )k 

with k = i(pi — p 2 ). In contrast to the conventional 
Skyrme ansatz, where the energy functional contains a 
Hartree- and a Fock-contribution, we here neglect the 
exchange (Fock-) term for the spin-orbit potential in the 
last line of Eq. (jlj). Otherwise, the operator P T would be 
equivalent to +1 for spin saturated systems. For the rest 
of the potential the Fock terms are included. The asso¬ 
ciated 11 parameters t. L , x t (i = 0... 3), Wo, x w and a of 
the Modified Skyrme Ansatz (MSkA) in Eq. ( 0 ) are de¬ 
termined by a fit to experimental data of finite spherical 
nuclei. The nuclear properties taken into consideration 
are the empirical binding energies and charge radii of the 
closed-shell nuclei 16 0, 40 Ca, 90 Zr, and 208 Pb. In order 
to take into account the variation in isospin we have also 
included Sn-isotopes 116 Sn, 124 Sn, and the doubly closed 
nucleus 132 Sn as well as one of the lead isotopes 214 Pb. 
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The resulting force and its parameters are presented in 
Table 

The spin-orbit term in the single-particle field derived 
^from this force has the form W T (r)(p x <x) with 

W T (r) = W.Vpr + ^Vpr^r, (2) 

where p T is the density for neutrons or protons (r = n or 
p) and W\ = Wo(l + x w )/2, W 2 = Wo/2. Conventional 
Skyrme calculations use a spin-orbit potential without 
isospin dependence and include the Fock term. This leads 
to x w = 1 and the relationship W 1 /W 2 = 2, which is 
nearly by a factor 2 different from the value 1.0005 ob¬ 
tained within the modified Skyrme Ansatz MSkA (see 
Table |). It is interesting to note that the fit leads to a 
value of x w very close to zero, which corresponds to the 
one without the Fock term. 

In order to study the spin-orbit term in the RMF the¬ 
ory we start from the standard Lagrangian density |^| 

C = i> (7 (id - - g p pf - eA) - m- g a a) ip 

+ \{8a) 2 - U(a) - + \m 2 u u 2 

+ ^ m y - ( 3 ) 

which contains nucleons ip with mass m. a-, to-, p- 
mesons, the electromagnetic field and nonlinear self¬ 
interactions U(a) of the cr-field, 

U(a) = ^mlcr 2 + ^g 2 <7 3 + ^g 3 a‘ l - (4) 

In a non-relativistic approximation of the correspond¬ 
ing Dirac equation [jlTj we obtain the following single 
particle spin-orbit term: 

W T (r) - +C 2 + C 2 p )Vp T 

m z m* z H 

+ + C 2 )V Pt ^ t , (5) 

m z m* z F 

with C 2 = ( rngi/rrii ) 2 for i = a,uj,p, which is similar 
in form to the spin-orbit field (i|), but which contains 
r-dependent parameters 

Wi = ^-^(C 2 + C 2 u + C 2 ) (6) 

ra z ra* z 

W2 = ^-^{C 2 a + C 2 -C 2 p ), ( 7 ) 

m z m z 

with m*( r) = m — g a v( r) The r-dependence drops out 
in the ratio W 1 /W 2 which is 1.13 for the parameter set 
NL1 H and 1.10 for the parameter set NL-SH Q. This 
is only slightly higher than the value 1.0005 obtained in 
the Modified Skyrme Ansatz by the fit to the empirical 
data (see Table |l[). The absolute size of the spin-orbit 
term turns out to be r-dependent, which stems from the 
r-dependence of the effective mass m*(r). It can be ap¬ 
proximated by m/m* «1 + C 2 /m 3 p( r). A more careful 


consideration would therefore require an explicitly den¬ 
sity dependent spin-orbit term. It has not been included 
in the present investigation. 

In Table || we show nuclear matter results obtained 
in the MSkA and compare it with the values from the 
conventional Skyrme force SkM* |l^]. The saturation 
density is obtained as po = 0.1531 fm -3 in MSkA. This 
is the same as the value obtained from an extensive fit 
of the mass formula FRDM ©• The binding energy 
per particle E/A is 16.006 MeV, which is close to that 
of other Skyrme forces. The compression modulus K = 
319 MeV is somewhat higher than that of the presently 
adopted Skyrme forces, but it lies within the error bars 
of the analysis based upon the breathing mode energies. 
The asymmetry energy J is close to the empirical value 
of 33 MeV. The effective mass m* is in good agreement 


with that of SkM*. 

In Table III we show binding energies and charge radii 
obtained in the MSkA for a number of spherical nuclei 
and compare them with those from SkM*. A comparison 
with the empirical values shows, that the binding ener¬ 
gies obtained with MSkA have improved over those of 
SkM*. The slightly reduced charge radii in MSkA seem 
to be connected with a slightly higher binding energy, 
which improves the results for the lighter Pb-isotopes. 
In general the charge radii are improved including that 
of 16 0. 


In Fig. 1 we show the isotope shifts of Pb nuclei for 
the Modified Skyrme Ansatz MSkA together with ex¬ 
perimental data and results of conventional Skyrme cal¬ 
culations. For MSkA we observe a clear kink at the 
double magic nucleus 208 Pb, whereas the conventional 
Skyrme force SkM* with an isospin independent spin- 
orbit term gives an almost straight line. For the lighter 
isotopes both theories give excellent agreement, on the 
heavier side MSkA comes closer to the experimental iso¬ 
tope shifts. It may be recalled, that the RMF theory 
||12| is successful in reproducing the full size of this kink. 
MSkA uses an isospin dependent, but density indepen¬ 
dent spin-orbit force. In contrast the spin-orbit term de¬ 
rived from the RMF theory (see Eq. |H|) is implicitly den¬ 
sity dependent through the density dependence of the 
effective mass m*(r). A density dependence of the spin- 
orbit term in Skyrme theory might improve the charge 
radii of the heavier Pb-isotopes. It has also been observed 
that a density dependent pairing force can possibly im¬ 
prove the situation in this context jli]]. This requires, 
however, further investigations. 

In Fig. 2 we present binding energies of Zr-nuclei about 
the neutron drip line as a function of the mass number. It 
is observed that in agreement with earlier investigations 
[fldl within the conventional Skyrme theory shell effects 
about the closed shell nucleus 122 Zr are weakened con¬ 
siderably. In contrast, RMF-calculations exhibit strong 
shell effects in the drip line region. It has been surmised 
in Ref. |l^| that this is caused by the differences in the 
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spin-orbit terms in the two approaches. It is gratifying 
to see that introduction of an isospin dependent (not an¬ 
tisymmetrized) spin-orbit term in Skyrme theory leads 
to stronger shell effects. This is also in agreement with 
the predictions of the FRDM |h| , where also an isospin 
dependent spin-orbit term is used. A more quantitative 
analysis shows that the ratio W\/Wi in FRDM is close 
to 1.06 for the Pb-nuclei and 1.09 for 122 Zr. 

Summarizing, we conclude that the isospin dependence 
of the spin-orbit term has an essential influence on the 
details of anomalous isotope shifts of Pb-nuclei. A new 
Modified Skyrme Ansatz has been proposed, and the 
isospin dependence of the spin-orbit strength has been 
determined. The magnitude of this isospin dependence 
x w is in agreement with the deductions from the relativis¬ 
tic mean-field theory. A reasonably good agreement with 
the experimental data on the binding energies and charge 
radii has been obtained. The kink in the isotope shifts 
of Pb-nuclei has been obtained in the modified Skyrme 
ansatz. However, the agreement with the empirical iso¬ 
tope shifts for heavy Pb-nuclei is not so good as in the 
RMF theory. This calls for further investigations includ¬ 
ing a density dependence of the spin-orbit term. 

Finally a remark about the isospin dependence of the 
spin-orbit term: our investigations lead to the interest¬ 
ing result, that the parameter x w in the Eq. (1) is close 
to zero. This leads practically to an isospin indepen¬ 
dent single-particle spin-orbit field in Eq. (2). This is 
in agreement with relativistic calculations, where the en¬ 
tire isospin-dependence of the spin-orbit field is caused 
by the parameter C p in Eq. (5), which is in fact rather 
small. In contrast the two-body spin-orbit potential in 
conventional Skyrme theory is isospin independent. The 
exchange term, however, causes a strong isospin depen¬ 
dence in the corresponding single-particle spin-orbit field. 
On the other hand in RMF theory the spin-orbit field has 
its origin in Lorentz covariance. There is no contribution 
from a two-body spin-orbit potential and an exchange 
term is therefore excluded. 
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FIG. 1. Isotopic shifts in the charge radii of Pb isotopes 
normalized to the nucleus 208 Pb as a function of the mass 
number A in the Modified Skyrme Ansatz (MSkA) as com¬ 
pared to the Skyrme force SkM* and the empirical data [fl9|. 


FIG. 2. Binding energy of Zr isotopes at neutron drip line 
in the MSkA and SkM* calculation. 


TABLE I. Parameters of the interaction in the modified 
Skyrme Ansatz (MSkA). 


to = -1200.074 (MeV fm 1 2 3 4 ) 

x 0 = 0.187 

ti = 396.302 (MeV fm 5 ) 

*1 = 0.018 

t 2 = -105.579 (MeV fm 5 ) 

X 2 = -0.059 

t 3 = 10631.527 (MeV fm 3+3a ) 

x 3 = 0.046 

Wo = 316.38 (MeV fm 5 ) 

x w = 0.0005 

a = 0.7557 



TABLE II. Nuclear matter properties obtained in the mod¬ 
ified Skyrme Ansatz (MSkA). 



MSkA 

SkM* 

Po 

0.1531 fm -3 

0.1603 fm -3 

C E/A )oo 

16.006 MeV 

15.776 MeV 

K 

319.4 MeV 

216.7 MeV 

III 

Po e 

100.2 MeV 

57.9 MeV 

J 

30.0 MeV 

30.0 MeV 

m* /m 

0.76 

0.79 
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TABLE III. The binding energies and charge radii ob¬ 
tained with the Modified Skyrme Ansatz (MSkA) in the 
Hartree-Fock approximation as compared with the normal 
Skyrme force SkM* and the empirical values. 


Nuclei 

Binding Energies (MeV) 

Charge Radii (fm) 

expt. 

MSkA 

SkM* 

expt. 

MSkA 

SkM* 

lb O 

-127.6 

-128.1 

-127.7 

2.730 

2.742 

2.811 

40 Ca 

-342.1 

-342.8 

-341.1 

3.450 

3.468 

3.518 

48 Ca 

-416.0 

-416.8 

-420.1 

3.500 

3.506 

3.537 

90 Zr 

-783.9 

-781.9 

-783.0 

4.270 

4.274 

4.296 

116 Sn 

-988.7 

-984.0 

-983.4 

4.626 

4.623 

4.619 

124 Sn 

-1050.0 

-1047.9 

-1049.0 

4.673 

4.677 

4.678 

132 Sn 

-1102.9 

-1106.3 

-1110.7 

- 

4.728 

4.727 

200 Pb 

-1576.4 

-1570.9 

-1568.4 

5.464 

5.465 

5.468 

202 Pb 

-1592.2 

-1588.3 

-1586.0 

5.473 

5.475 

5.478 

204 Pb 

-1607.5 

-1605.4 

-1603.4 

5.483 

5.486 

5.489 

206 Pb 

-1622.3 

-1622.3 

-1620.3 

5.492 

5.496 

5.501 

208 Pb 

-1636.5 

-1637.8 

-1636.4 

5.503 

5.506 

5.510 

210 Pb 

-1645.6 

-1645.8 

-1645.6 

5.522 

5.522 

5.520 

2i2p b 

-1654.5 

-1653.8 

-1654.5 

5.540 

5.537 

5.531 

214 Pb 

-1663.3 

-1661.6 

-1663.1 

5.558 

5.552 

5.541 
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